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CT -1 Abstract 

We study a porous medium equation with nonlocal diffusion effects given 
43 by an inverse fractional Laplacian operator. More precisely, 

p u t = V ■ (uV(-A)- s u), 0<s<l. 



The problem is posed in {x G R N ,t € M} with nonnegative initial data u(x,0) 
that are integrable and decay at infinity. A previous paper has established 



the existence of mass-preserving, nonnegative weak solutions satisfying energy 
estimates and finite propagation. Here we establish the boundedness and C a 
\q regularity of such weak solutions 
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1 Introduction 



This paper is devoted to study the regularity properties of weak solutions of a model 
of porous medium equation that includes nonlocal effects through an integral relation 
of pressure to density. This allows to account for long-range effects. 

Let us recall the typical derivation of the porous medium equation, cf. [TJ [19]. 
We consider a gas propagating in a homogeneous porous medium; its dynamics is 
described by first assuming conservation of mass 

d t u + V ■ (vu) = , 

where u(x,t) > denotes the density of the gas and v(x,t) is the (locally averaged) 
velocity. We then postulate that the motion proceeds according to Darcy's law so 
that v = —Vp, where the velocity potential is interpreted as a pressure. Finally, some 
barotropic state law for gases implies that p is a monotone function of u, p — f(u). 
In this way we get the equation 

dtu = V-(«V/(it)). 

The simplest case (called isothermal) is p = u and in that case we arrive at dtu = 
cAu 2 , which appears in a different context as a model in groundwater infiltration, 
Boussinesq's equation [21 H] . 

The novelty of our present model consists is relating p to u through a linear integral 
operator that makes a kind of average of the space distribution u(-,t), 

(1.1) p(x,t) — Cu(x,t), Cu(x,t) := J L(x — y) u(y,t) dy 

More in particular, the positive kernel L is locally integrable and decays slowly at 
infinity to represent "long-range " interactions. To be specific, we will work in WL N , 
we will take L(x) = c\x\~ N+2s , which is equivalent to saying that p is given as an 
inverse fractional Laplacian, i.e., p — (—A)~ s u, and we consider < s < 1. 

In a previous paper [11 J we have introduced this model and proved existence of weak 
solutions for the Cauchy problem 

(1.2) u t = V-{uVp), p = Cu = {-A)- S u, < s < 1 , 
posed for x G Mr, N > 1, and t > 0, with initial conditions 

(1.3) u(x,0) =u (x), xeR N , 
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where uq is a nonnegative and integrable function in M. N decaying as |x| — >• oo. 

Let us point out that equations of the more general form u t = V ■ (cr(u)V Cu) have 
appeared recently in a number of applications in particle physics. Thus, Giacomin 
and Lebowitz consider in [13] a lattice gas with general short-range interactions and a 
Kac potential J 7 (V) of range 7 -1 , 7 > . Scaling spacelike with 7 _1 and timelike with 
7~ 2 , and passing to the limit 7 — )• 0, the macroscopic density profile p(r,t) satisfies 
the equation 



(1-4) t 



, ,„6F(p) 



Here F(p) = J f s (p(r))dr — (1/2) J J J{r — r')p(r)p(r') drdr' , where f s (p) is the 
(strictly convex) free energy density of the reference system, and cr s (p) is the mo- 
bility of the system with only short-range interactions. See also [14] and the review 
paper [JS|. The model is used to study phase segregation in 



Further motivations for model (1.3) can be found in [IT] and [21] , which contain 
references to applications in dislocation dynamics and in superconductivity, as well 
as current mathematical progress. 

Mathematical results. Paper [H] contains the proof of existence of a weak solu- 



tion of Problem ( 1.2 )-( 1.3 ) when uq is a bounded function and has exponential decay 
at infinity. Besides, a number of basic properties are proved, like energy estimates, 
bounds in the L p spaces, and the property of finite propagation that says that com- 
pactly supported data produce solutions whose support is compact in space for every 
positive time. However, the question of uniqueness of weak solutions is a pending 
open problem (in more than one space dimension). Comparison theorems, a crucial 
tool in parabolic equations, are only available under special circumstances (i.e., for 
so-called true super- or sub-solutions). The asymptotic behaviour of the solutions as 
t — y 00 has been studied by two of the authors in [TJ] using obstacle problems and 
entropy estimates. 

The regularity theory that we develop in the present paper is as follows: 

a) If u has initial data in L 1 logL, then it becomes instantaneously bounded; 

b) If u has initial data in L 1 , then it falls into the previous case 

c) Bounded solutions are continuous with a modulus of continuity. 

After some preliminaries, Section [2j and the needed theory on bilinear forms con- 
tained in Section [3j the boundedness results are stated as Theorem 4T and proved 



in Section |4j The proof of the C a regularity result, Theorem 5.1 takes up Sections 
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[5jto 12 It says that for < s < 1, with s ^ 1/2, bounded solutions u > are C a 
continuous in (x,t) with some universal exponent a G (0,1) that depends on N,s. 
The proof of this result in the range < s < 1/2 uses a number of techniques that 
are becoming classical in the study of regularity of nonlocal diffusion problems, but 
it is complicated since we must take into account both the nonlinearity and the pos- 



sible degeneracy. Section 12 covers the more difficult range s > 1/2. The regularity 
result in that case uses transport ideas in the form of a geometrical transformation 
to absorb the uncontrolled growth of one of the integrals that appear in the iterated 
energy estimates. The case s = 1/2 has new difficulties and will be treated separately. 



As a consequence of these results, in Section [13] we complete the existence theory by 
constructing a continuous weak solution for any initial data u Q e L 1 (IR Ar ), u > 0. 



Notations. We will refer to Equation (1.2) as the FPME (with F for fractional). We 



will use the notation (— A) s with < s < 1 for the fractional powers of the Laplace 
operator defined on the Schwartz class of functions in M. N by Fourier transform and 
extended in a natural way to functions in the Sobolev space H 2s (M. N ). Technical 
reasons imply that in one space dimension the restriction s < 1/2 will be observed. 
The inverse operator is denoted by C s = (— A)~ s and can be realized by convolution 



(1.5) C s u = L s -ku, L s (x) = c(N, s)\x\ 



-N+2s 



1/2 

C s is a positive self-adjoint operator. We will write V, s = C s which has kernel L s / 2 . 
The subscript s will be omitted when s is fixed and known. For functions that depend 
on x and t, convolution is applied for every fixed t with respect to the space variables. 
We then use the abbreviated notation u(t) = u(-,t). 

For a measurable u > and for k > we denote by = (u — k) + = max{« — k, 0}, 
and = min{w — k,0} in such a way that > > u^, the supports of u\ and 
agree only on points where u = 0, and also u = k + + . We will use similar 
notations: u+ = (u— y?) + , u~ = (u — ip)~ when ip is a function and not just a constant, 
and then we may split u as follows: u = (p + u+ + u~ . 



2 Preliminaries. Existence and basic estimates 



Definition. We say that u is a weak solution of Problem (1.2)-(L3j) in Qt 



(0,T) with initial data u e L 1 ^) if u e L\Q T ), C{u) e L} oc {U,T : W^(R N )), 



X 
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and uVC(u) G L 1 {Qt), and if the identity 
(2.1) u (r] t - V£(u) ■ Vt]) dxdt + I u (x) r)(x, 0) dx = 



holds for all continuously differentiable test functions rj in Qt that are compactly 
supported in the space variable and vanish near t = T. 

The following results have been proved in 

Theorem 2.1 Let u G L°°(R N ), u > 0, and such that 
(2.2) u (x) < Ae~ alxl for some A, a > . 



Then there exists a weak solution u of Equation (1.2) with initial data Uq. Besides, 
u G L°°(0,oo : L 1 (R N )), u G L°°(Q), VU{u) G L 2 (Q). Moreover, for all t > we 
have conservation of mass: 

(2.3) / u(x,t)dx= / uq(x) dx , 

Jr n Jr n 

as well as the L°° bound: ||w(t)||oo < ||mo||oo- The constructed solution decays expo- 
nentially as \x\ — > oo. The first energy inequality holds in the form 

(2.4) f [ \VHu\ 2 dxdt + [ u{t) log(«(t)) dx < [ u log{u )dx, 
Jo Jr n Jr n Jr n 

while the second says that for all < t\ < t 2 < oo 

(2.5) f [ u\VCu\ 2 dxdt + \ [ \Hu{t 2 )\ 2 dx <]- [ [H^uit^ 2 dx . 

It, JwN 2 ./iffAT 2 ImN 



Other properties of the constructed solutions. Here are some of the most useful 

• Translation invariance. The equation is invariant under translations in space and 
time, and this property reflects on the set of weak solutions. 

• Scaling. Moreover, the equation is invariant under a subgroup of the group of dila- 
tions in (u,x,t), and this implies a scaling property for the set of solutions. Namely, 
if u(x,t) is a weak solution as described in the existence theorem, with initial data 
uo(x), and A, B, C are positive constants, then u(x, t) = Au(Bx, Ct) is again a weak 
solution on the condition that A = CB~ 2+2s . It has initial data uq{x) = Auo(Bx). 



5 



• Conservation of sign: uq > implies that u(t) > for all times. 

• LP estimates. The LP norm of the solutions, 1 < p < oo, does not increase in time. 

• Finite propagation: Compactly supported initial data u (x) give rise to solutions 
u(x, t) that have the same property for all positive times, i.e., the support of «(•, t) is 
contained in a finite ball Br^(0) for any t > 0. 

• A standard comparison result for parabolic equations does not work in general. 
This is one of the main technical difficulties in the study of this equation. In fact, 
special situations are found in [UJ where some comparison holds by using so-called 
true super- and sub-solutions. 

Energy solutions. The constructed solutions are limits of smooth functions for 
which the energy inequalities are justified. In the sequel we will need this fact and also 
similar integrations by parts involved in the new energy inequalities. In particular, 
we want the weak solution to satisfy the identities 



u (f] t — B r (u, 1])) dxdt + / uq(x) r](x, 0) dx = 



(2.6) 



where B r is the bilinear form that will be defined in the next section, r 



s and 



1] E L 2 (0,T : LL r (M. N )), r\ bounded, r\ t E L 2 {Q T ). This class of solutions can be called 



weak energy solutions. Below (see Formulas (4.4), (4.6)), we will need a version a 
this definition what consists in using 7] = f(u) and integrating in time to get 



(2.7) 



F(u(t))dx 



tl 



+ It I v uvcu dxdt = °' 



with / is smooth and bounded, F(s) = f s f(s) ds, and < t± < t% < T. Note that 
we do not need to assume regularity for r] t . The constructed bounded solutions are 
energy solutions in this sense. We will also use t] = f(u/(p), where / as before and 



tp(x) is a smooth positive function that does not vanish, see Formula (6.2) and later. 



3 Bilinear forms 



Before proceeding with the study of the boundedness and regularity properties, we 
need some results on fractional operators. The bilinear form associated to the space 
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H r (R N ) = W r > 2 (R N ), < r < 1, is 

(3.1) B r (v,w) = C N J J (v(x) - v(y)) . _ \ N+2r {w{x) - w(y)) dxdy 

It is easily seen by Fourier transform that this is well defined for functions in H r (R N ). 
We will omit the subindex and write B instead of B r when the context is clear. 

Corollary 3.1 (a) If v is a monotone function of w, i.e., if v = G(w) with G' > 
then 

B(v,w) > 0. 

(b) If G'(s) < C for some constant C > 0, then 

B(G(w),w) < CB{w,w). 

Proposition 3.2 For every u,v G H 1 (R N ) = W 1 ' 2 (R N ) we have 

(3.2) B r (v,w) =\\ Vv(x)-. l—-Vw(y)dxdy. 



\x-y 



Proof. Prove first for functions and use a smoothing and truncation of the kernel. 
Then pass to the limit. □ 

Remark. Since the weak formulation of the FPME leads to an expression of this 
latter form with kernel K = c\x — y\~ N+2s , we will put below r = 1 — s. Actually, 
all that we will use in Sections [4] and later, in accordance with Proposition 3.2 



is 



two kernels L and K such that L, K > and AL = K, as well as the bilinear forms 
associated to the pairings 

V(-) ■ L ■ V(-) < — > difference ■ K ■ difference, 



which is a short way of writing the equivalence of formulas (3.1), (3.2). In later 
calculations we will also use the following positivity properties of the integration of 
these kernels applied to truncations of functions. 

Lemma 3.1 Let u e H r (R N ) and u\ G H r (R N ). Then, 
(3.3) B r (u+,u)>B r (ut,u+). 
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Proof. If K r (x,y) — \x — y\ ( N+2r ) is the kernel of B r , we have 

B r (u+, u) - B r (u+, u+) = B r (u+, u - u+) = B r (u+, k + ul) = 
[u+ (x) - u+ (y) ) K r (x, y) (u^ (x) - (y) ) dxdy . 



Now, given the fact that u + (x)u (x) = a.e. and symmetry in x, y, the last integral 
equals 

-2 J u%(x)u^(y)K(x,y)dxdy > 0. □ 
We will also need the following embedding inequality. 
Lemma 3.2 For every u G L\R N ) n H r (R N ) we have 

(3.4) J u q dx < C\\u\\ e Ll \\u\\ 2 Hr , 

where 9 = 2r/N , q = 2 + 9 = 2 + (2r/N), and C = C(N, r) > 0. Moreover, for every 
u E L 2 (R N ) n H r (R N ) we have 

(3.5) J u qi dx < C\\u\\ e L2 \\u\\ 2 Hr , 

where 6 = 4r/N, q l = 2 + (Ar/N), and C = C(N, r) > as before. 
Sketch of the proof. We use the Sobolev inequality that says that 

(3.6) (^J u p dx^j ' < C\\u\\ 2 Hr 

for some p > 2 depending on r G (0, 1) and N. C depends also on r and N. Actually, 
p = 2N/(N - 2r), when iV > 2 or when TV = 1 if, in addition, < r < 1/2. We want 
to control for some q > 1 using that 

u q dx < ( I udx ) ( I u p dx 



where q = 9 + p(l — 6). We will take the values 

p p 



The proof when u G L 2 is quite similar. 

The case N = 1, with l/2<r<lis easy. Take < r' < 1/2 and observe that 
L 1 fl H r C H r as a continuous embedding. Now we can use (3.6), with r' replacing 



r, and get (3.4) with, perhaps, some different values of q and 9. □ 



Remark 3.3 We also recall that for every convex function $ the quantity J $(«(£)) dx 
is non-increasing in time in the FPME evolution. We note that for such a function 
$ we have 



$(u(x,t)) dx 



Ti 



+ I B(G(u),u)dt < 0, 

T J T 

where G(u) is the primitive of§"(u)u. The bilinear form B is as above and 



B(G{u),u) > cB{u,u) 



\u 



if G' is strictly positive. We will use this in the case of the truncations in the form 
with ul = {u — k) + and then G = for u < k. 



4 Boundedness of solutions 

This section is devoted to proving the main boundedness result. Here < s < 1. 



Theorem 4.1 Let u be a weak energy solution of Problem ( 1 . 2 ) — ( 1 . 3 ) with uq £ 
L 1 (IR iV ) nL°°(]R Ar ) ; as constructed in Section 2. Then, there exists a positive constant 
C such that for every t > 

(4.1) sup \u(x,t)\ < cr a \\u Wl 1(Rn) 

with precise exponents a = N / (N + 2 — 2s), 7 = (2 — 2s) /(N + 2 — 2s). The constant 
C depends only on N > 1 and s G (0, 1). 

In dimension N — 1 this is proved in [3]. Our proof applies to all iV > 1 and is 
divided into three subsections. Note that this estimate and conservation of mass 
imply a decay for all intermediate norms L p with 1 < p < 00: 

(4-2) \H;t)\\ p <C p t-^\\u \\^ {Rn) , 
where a p = a{p — l)/p and j p = (1 + j(p — 
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4.1 Better integrability properties for solutions 

We start with the following partial result. 

Lemma 4.1 ("From L l to LlogL") Let u > be a weak energy solution of the 
FPME. If the initial data are integrable, then u(-,t) e LlogL for all positive t > 
and for all small t <t Q we have 



(4.3) J u(t) log(l + u(t)) dx + - J sB(u(s)t,u(s)t)ds < C | log(t)|||«o||i , 

where Co depends only on N,s, and we write u(t)f = (u(t) — 1)+. The time to is 
estimated as to = inf{l, C||uo||^i } for some constants C, i? > 0. 

Proof. We use as test function rj = £log(l + u) on the weak form of the equation. 
After observing that h(u) = (u + 1) log(l +u) > satisfies log(l +u) — h'{u) — 1 and 
integrating in M. N x (0, r), we get the identity 
(4.4) 

r h(u(r)) dx + / tB(g(u),u) dt = / / h(u(t))dxdt + r u(r) dx — // udxdt, 



where B = B r with r = 1 — s, as already explained, and g'iu) = u/(l + u) with 
g(0) = 0. The last two terms in the display disappear by mass conservation. Note 
next that g'(u) > 1/2 for u > 1, hence ^(u) > (1/2) (w — 1) + , so that, writing 
uf = (u — 1) + we have by the already monotonicity properties of B: 

B{g{u),u) > (l/2)B(u+,u) > (l/2)B(u+ 

With this we arrive at 

(4.5) r [ h(u(r))dx+- [ tB(uf,uf)dt < [ [ h(u{t))dxdt . 



■2 

Recalling the definition of h and using log(l + u) < u, we get 

sup t I u{t) log(l+w(t)) dx+- I tB{u[ ,u[) dt < f f u\og{u+l) dxdt- 
{o<t<r} J 2 y Jo J 



+r / u dx. 



We still have a "bad term" in the right-hand side containing J wlog(l + u) dx, and 
it is tackled as follows: note that the left-hand side controls J tdt J u q dx if q is as 
in Lemma 3^ take such a q > 2. Then, for any M > 2 we have 

I u log(l + u) dx < log(M + 1) ! udx+ l ° g ^ + _^ I u q dx , 

J Ju<M M q Ju>M 
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where we have used the fact that u log(l + u)/u q is decreasing in u for u > M. Next, 
there is a constant C(N, q) > such that u q < C(u — l) q + for this range of u. Choose 
now M = t~ a and then (3 < a(q — 1) — e (taking care that (3 > 1). With all this, we 
get for small r 

f I ' u log(l + u) dxdt <d f dt\ log(t) | ( sup f u)+C 2 f t p dt{ f (u - 1) 9 + dx) 
o J Jo t J Jo J 

<dT\]oe{T)\ [ U dx + C 2 r^\\u(t)t\\U\u+(t)\\ 2 H rdt. 



Since (3 > 1 the last term is controlled by the 0-energy term in the left-hand side of 
(4.5) for small r. In particular, we choose < r < 1 and t 13 ^ 1 < (2<Z7 2 1 1 iio 1 1 rO -1 - The 
other term is a multiple of f u , hence bounded. We get 

sup t / u(t) log(l + u(t)) dx + / tB(u± , ) dt < C 3 r| log(r)| ||wo||i- 

{0<t<r} 7 Jo 

Putting t = t we get the result. □ 

Lemma 4.2 ("From LlogL to L 2 ") Initial data in the space LlogL imply that 
u(-,t) G L 2 for all positive t > and the bound on the L 2 norm of {u — 1) + de- 
pends only on t, s,N, ||wo||i; an d \\( u o ~ l) + IUio g L- 

Proof. We define v(x, t) = u(x, t) V 1 so that v = 1 + uf > 1 and v t = u t x{u > 1). 
Recall the notation uf = (u — k) + . Then, 

d 
dt 



v \og(v) — v] dx — J \og(v ) v t — J log(t> ) u t . 



Using the weak form of the equation with rj = log v as test function, we get 

*2 r t 2 



(4.6) / [v\og(v) ~v}(t)dx 



V[log(u)] uVCudxdt = 0. 



We work out the last term for fixed time and observe that, since u — v for u > 1, we 
have 



/ V[log(v)] uVCudx = I Vv- VCudx = / Vuf ■ VCudx = B r (uf, 

Ju>\ J J 



U 
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Using again the monotonicity of B, see Corollary 3A and Remark 3^3 and putting 
h(v) = v log(f ) — v, we get 



h(v(t 2 ))dx+ / B r (uf,uf)dt< / h(v(ti))dx. 



(4.7) 



Note that h(v) is convex for v > 1 and the right-hand side is bounded by a combination 
of ||w ||i and ||(w - l) + |Ui og L- Hence, 



/ 2 B r {ut,ut)dt < C. 
Jo 



~ \\u 



1 K i ")\\HTl 



with r = 1 — s. Use Lemma 3.2 



to 



(4.8) 

Recall finally that B r (u^,u^ 
conclude that uf(t) G L q (M. N ) for some q > 2. More quantitatively, this together 
with (4.8) and Remark 3.3| give the estimate 

supt\\ut(t)\\ q Lg <C. 



(4.9) 



i>0 



with C depending as in the statement of the lemma. Interpolation with L 1 gives the 
result. □ 



4.2 Boundedness 

With the preceding results, we may assume that uq G L 1 (IR iV ) R L 2 (IR Ar ) after some 
displacement of the time origin. Then we can follow the De Giorgi approach (as 
outlined for instance in [9], [ID])- We consider the truncations 

uj(x,t) = (w-M(2-2- J ))+ 

The value of constant M > will be conveniently chosen later. Actually, we may 
assume J dx very small by selecting M large. 

Claim. The following energy inequality holds for all < t\ < t^- 

(4.10) J uf{t 2 ) 2 dx + 2M(2 - 2~ j ) ' B{u+ , w+) dt < j u+ih) 2 dx . 

To see this, we use the definition of weak solution for our FPME with rj = as a 
test function. Then, for t\ < t 2 , we have 

J u t u+dxdt=^ J u+(t 2 ) 2 dx-^ J uj(t 2 ) 2 dx. 
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For the RHS we observe that u = + M(2 — 2 ■?) whenever u > M(2 — 2 J ). Hence, 

1 



tt>Af(2-2-J) 



J 2 



V(w+) 2 - V£wdx + M(2-2^') / Vut-V£uda; 



£((w+) 2 , u) + M(2 - B(u+ u) > M(2 - 2" J ') B(uj, 



In the last inequality we have used both Corollary 3.1 and Lemma 3.1 This gives 

Q . □ 

We now fix to > 0. We want to prove that the solution it is bounded for all times 
t > t . As in [TU], let us define the total energy for the truncated function u- as 

//*oo 
{u+) 2 {t)dx + 2M I B{u\,u\)dt, 
JTj 

where Tj = t (l — 2~ J ). From (4.10), taking arbitrary values t 2 = t > Tj and 
ti = t' G [Tj_i,Tj] we have 



(4.11) 



Aj < inf 



(«+) 2 (0^. 



Observe now that u^{x) > implies ■u^_ 1 (x) = u^(x) + M2 5 > M2 ■ 7 . Therefore, 
for every p > 2 we have (keeping the time fixed) 



da; = / («+_! - M2^) 2 • X{u+>0} dx 



p-2 



2-J'M 



MJ 



p-2 



If p > 2 is the exponent corresponding to Sobolev's embedding theorem, we deduce 

2 j \ P-2 



from (J6U21 that 



< C JV (2VM) p - 2 inf t , e[T3 ._ 1)T3 . ] x ){t'f p/2 
Taking averages in t' we arrive to the inequality 



Aj < C N (2?jM) 



p-2 



Tj - Tj-i 



T, 



Biu+^u+^mdt 



3-1 



p/2 



13 



< c 



N 



2i(p- 2 )2-?p/ 2 
MP- 2 (Mt )p/ 2 



M 



TP/2 



2-1 



This leads to a recurrence relation of the form 



Aj-x 

Mtn 



l+S 



with 5 



— 1 > 0, that implies Aoo = if Ao/Mto is small. This determines 



the correct value of M to choose. The conclusion is then that u(x, t) < 2M for all 
t>t . □ 



4.3 End of proof of the theorem 

The preceding subsections have established the result for any fixed t > 0, and we 
know that 

(4.12) u(x,t) <C(N,s,\\u \\i,t) 

but we do not know the dependence of C on its arguments in a precise way. We need 
to prove that this dependence takes the form (4.1 ). This is just a consequence of the 
scaling group that allows to pass from a solution u(x, t) to the rescaled solution 

(4.13) u(x,t) = Au{Lx,Tt) 

on the condition that A = TL~ 2+2s . On the other hand, we want to reduce u to unit 
mass, f u(x,t) dx = 1, and this means A = L N /||«o||i- All together this gives (with 
K||i = M) 

L = {MTf A = m^-^T"^, 
where (3 = 1/(N + 2 — 2s). We now apply the boundedness result to u at t = 1, i. e., 

u(x, 1) < C(N, s) VxGl N 

Going back to the u, we have 

su P m(x,T) = A sup«(a;,l) = M 2{1 - s)p T~ Np C{N, s) 

X X 

which gives the desired result upon replacing T by t. For other instances of this 
scaling argument cf. (20]. □ 
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5 Holder regularity. Main result and basic lemmas 



Once the question of boundedness is settled, we proceed with the local regularity of 
solutions. This is the main result. 

Theorem 5.1 Let u > be a bounded weak energy solution defined in a space-time 
strip S = R N x [Ti, T 2 ] C R N+1 . Let s e (0, 1), s =fi 1/2. Then u is C a continuous 
in the interior of Q for some exponent a(N, s) G (0, 1) and a constant that depends 
also on the dimensions of the subdomain Q <s S and the bounds on u in L°°(S) and 
L°°((Ti,T 2 ) : L 1 (R N )). 

Strategy. Since the equation is space- and time-invariant we may assume that 7\ < 
T2 = 0, and then we may study the regularity around x = and t = 0. 

The main ideas are two: on the one hand, we will prove some basic De Giorgi-type 
oscillation lemmas that say that the oscillation of the solution u decays when we 
restrict a basic domain, say, the cylinder T4 = [—4, 0] x -84(0), into a smaller cylinder 
like Ti = [—1, 0] x £>i(0). The second ingredient is the scaling property of the equation 
that allows to renormalize the solution through the transformation 

(5.1) u(x,t) = Au(Bx,Ct) 

with C, B > free parameters, and A = CB~ 2+2s . The way of attacking the problem 
is through the iterated application of the lemmas. At the end of every step we 
renormalize the solution defined in Ti into a rescaled solution defined in and we 
start a new application of the oscillation lemmas. In this way, we will show that 
the oscillation of the solution u decays dyadically in a family of space-time cylinders 
shrinking dyadically to a point. 

The needed lemmas have a simpler expression for < s < 1/2 where the diffusion is 
more similar to the standard porous medium case. For 1/2 < s < 1 convection effects 
appear that make some integrals diverge, and this makes the analysis more difficult, 
needing new techniques. The detailed study of how to proceed in the case s > 1/2 



are contained in Section 12 Until then we assume that s < 1/2 



5.1 The oscillation reduction lemmas 

These technical results need only be proved for bounded nonnegative weak solutions 
defined in a strip Sr = [— R, 0] x WL N . We denote by Tr the parabolic cylinder 
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[— R, 0] x Br(0). One of the lemmas controls the decrease of the supremum of the 
solution once we restrict the size of the parabolic neighborhood of (0,0), the other 
one implies that under suitable assumptions the solution separates from zero. A third 
one improves the first result so as to obtain a real alternative between going a bit 
down and a bit up. which leads to the proof of regularity. Here is the first basic 
lemma. 

Lemma 5.1 Let < s < 1/2. Given u, G (0, 1/2) andeo small enough {in particular ,eo < 
1 — 2s) , there exists 5 > (depending possibly on fi, e , s, and N) such that if we as- 
sume that 

(i) the solution u is bounded above in the strip S4 = W" 1 x [—4, 0] by 

(5.2) V(x) = 1 + (|x| £ -2)+, 0<e<e , 

and (ii) u is mostly below the level 1/2 in T 4 = £4(0) x [—4,0] in the sense that 

(5.3) \{u > 1/2} nr 4 | < 5|r 4 | , 

then we can lower the upper bound inside a smaller cylinder in the following quanti- 
tative form: u\ T < 1 — /i. 

We summarize the result by saying that "being mostly below 1/2 in space-time 
measure pulls down the supremum in a smaller nested cylinder". Note that for this 
lemma 8 can be chosen as a non-increasing function of /i with 5(1/2+) = 0. Also, 
e can be chosen as small as we want by sacrificing the gain in oscillation. We also 
remark that the size of the cylinders can be changed, though this affects the values 
of 5 if the new sizes do not conform with the parabolic scaling. Finally, the levels 
u—1/2 and u — 1 are taken by convenience, any pair of levels < Mi < M 2 will do, 
though in principle the value of 5 will change. 

A similar result applies from below but the proof is different since the equation is 
degenerate at u — 0. The idea is that if u is very often far from zero in T4 then in 
a smaller, suitably nested cylinder u stays uniformly away from zero. The technical 
version explains how "most of the space-time above 1/2, pulls up away from zero". 

Lemma 5.2 Under the same assumptions, given u, G (0,1/2) there exists 5 > 
(depending possibly on /U ,e 0! s, and N) such that if u satisfies 

(5.4) \{u> i/2}nr 4 | > (i-<s)|r 4 |, 

then u\ T > fiQ. 
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Again, 5 is a non-increasing function of /lq. Let So = 5(1/4), that is, when /lq = 1/4. 
A more elaborate version of this lemma will be needed in some cases of the iteration. 

We complement these two lemmas with a lemma that replaces the sentence "most 
of the space-time" of Lemma |5.1| by "in some set of positive measure" . 



Lemma 5.3 ("Some of the space-time below, pulls down") Assume as before that 
< s < 1/2 and u is trapped between and $ in S4. Besides, assume now that 

(5.5) \{u< i/2}nr 4 | ><y |r 4 |, 

with 5q defined as above. Then u\ r < 1 — /J , for some //(do). 

Notice that this second lemma applies only in one direction, reducing the oscillation 
from above. As in the classical porous media, we cannot expect this lemma to hold 
in the "pulling-up" case, due to the property of finite propagation (existence of so- 
lutions with compact support), a consequence of the degeneration of the equation. 
Nevertheless, this one-sided improvement will be enough to prove that the oscillation 
decays dyadically as follows: 



The iterated use of the Lemma 5.3 from above after rescaling at every step, as long 
as it is possible, reduces the oscillation of u from above and we start iterating and 
renormalizing to get U2,us, ■ ■ ■ . We note though that, as we do that, a renormalized 
solutions develop a "tail" in the sense that the functions Uk start to grow at infinity by 
effect of the scaling (both in vertical and horizontal directions). This is the reason for 
the form of the upper barrier that we use, which has an e-tail. Indeed, after k steps, u 
will be bounded by (1 — /i)~ k outside the B~ k dilation using (5.1). This first difficulty 



can be dealt by playing the integrability of the kernels L or K at infinity against a 
slow power growth in u. Indeed, by sacrificing the modulus of Holder continuity we 
may assume that the gain is very tiny (i.e., /i very small). Then, the build up in 
u(y,t) as y tends to infinity will be very slow (like \y\ £ ), being absorbed by K. If 



Lemma 5.3 never fails along the iteration, then we are at a point where u = and a 
Holder exponent is also found. 

On the other hand, if the process breaks down, then the first time that Lemma 



5.3| fails, it puts us into the hypotheses of Lemma 5.2 and that pulls u away from 



zero by a fixed amount /j,q (the need for an alternative is what makes Lemma 5.1 



insufficient). Then the operator becomes non- degenerate in the subsequent iterations 



and a counterpart of Lemma 5.3 applies also from below (upwards) since u will always 



be bounded between /i and 1. From there on the gain on the oscillation of u may 
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come from above or below in the dyadic way we have shown. The details of the 
iteration are given in Section 11 after the lengthy and delicate proof of the Lemmas 
is completed. We recall that all this will be done below for < s < 1/2. 



6 Lowering the maximum. Proof of Lemma |5.1 



We start here the technical work. The basic idea in the proof of the result is a 
particular kind of "localized energy inequalities" that will be iterated in the style of 
Giorgi to obtain the reduction on the maximum in a smaller domain. Localization 
is obtained by using a suitable sequence of cutoff functions. In order to deduce the 
necessary energy inequalities we use integration by parts formulas and analysis of the 



kernels. A main role is played by the bilinear form B(v,w) as defined in (3.2) with 



kernel K(x) = c\x\ ( N+2r ) and r — 1 — s. Moreover, we put L{x) = C\\x\ N+2s so that 



-AL = K. We will repeatedly use the equivalent form justified by Proposition 3.2 



Lemma 6.1 Let u be smooth and compactly supported, let v have small growth at 
infinity. Then, 

(6.1) B(u,v) = J j \u{x) — u{y))K{x — y){v{y) — v{x)) dxdy. 

We take a weak energy solution defined in a strip Q = (-T,0) x R N in the sense of 
Section |2j We justify the computations by recalling that u can be approximated by 
smooth positive solutions of similar problems as done in [TT] . 

[6[l. An energy formula. We consider a sequence of cutoffs tp(x) that have the 
form of perturbations of the level u = 1 within a region containing the unit ball -Bi(O), 
and an "outer wing" rising up above the 1-level for larger values of \x\. An explicit 
choice suiting our purposes will be done below. We only need to know at this stage 
that the cutoff function tp is smooth, lies above 1/2 everywhere, and also that u < tp 
for all \x\ > 3 for all times —4 < t < 0. 

We use the function rj = \og((u/ip) V 1) = log(g) as a test function in the weak form 
of the equation. Note that 

g : = %l = l+ ( "- y)+ = l4, 
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where ut = (u — p) + according to the adopted notation. Note that g > 1 and g > 1 
iff u > <p. According to our assumptions, w+ and g — 1 have compact support in the 
ball of radius 3. We will often split u as follows 

u = u+ + (f + (u - (f)~ 

where we write (u — ip)~ = (u — <p) A = u~. Notice that with this notation we have 
u~ < 0. After applying the weak formulation of the equation with r\ as above, we get 
on the LHS for T x < t < T 2 < 0: 



(6.2) 



T 2 
Ti 



rjd t udxdt — j J log ^— V 1^ ipd t (u/(p) dxdt — 









T 2 f 


") 


log 


dx 










Ti J 



v(g logs- -g)dx 



r 2 



We will need an estimate of this quantity: after adding 1 to the last integrand we get 
the expression H(g) := g \ogg + 1 — g for which we have the estimate H(g) ~ (g— l) 2 
for for 1 < g < 2, in the sense that 



(6.3) 



1 fu 
2 



^ < H(g) < 
V J \V 



u2 



Let us now calculate the right-hand side of the expression in the weak formulation 
of the equation. We have 

dt \og(g(x)) div [u(x) VL(x — y) (u(y) — u(x))] dxdy 

( 6 - 4 ) r f Vg(x) 

dt / — j—— u{x)V x L{x — y) (u(y) — u(x)) dxdy — I + II , 
Ju>ip 9\ x ) 



u ul 

g= - VI = 1 + 



where we pass from the first line to the second integrating by parts. Recalling that 

Vg = V(u+/p)x({u>p}), 
the first part of the splitting is: 

1 = ~J dt J J u+{x)(-/\L)(x -y)[u(y) -u{x)}dxdy 
dt J J u+{x) K(x - y)[u{x) - u{y)} dxdy . 
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After symmetrizing, we get I = — | J B(u+, u) dt, where 

£(uj, u ) ■= J J fa J 0*0 ~ K(y)) K ( x ~ y)( u ( x ) - u (y)) dxd v- 



On the other hand 
II = 



dt J J u+(x) ^^ V x L(x - y)[u(y) - u(x)]dxdy := J Q(u+,u) dt. 

• In order to separate the good and bad components of both B and Q, we use the 
decomposition u = m+ + cp + u~ . We get 

»(«+«) = u+) + B(«+ ^) + u-). 
and a similar expression for Q. We now make some observations: 

(i) B(u+,u+) is a positive quadratic form. We will pass the corresponding part of I 
to the LHS as a term with positive sign and thus complete the energy expression in 
the energy inequality that we want to derive. 

(ii) B(u^,u~) has also the correct nonnegative sign because of these facts: u+ and 
u~ have opposite signs and disjoint supports, and K > 0. We could drop this term 
in a first calculation, but we will keep it and use it to control some of the bad terms 
in Q. 

Summing up, we have up to now the basic identity for T% < T 2 < 0: 



(6.5) < 



^(glogg + l- g)\ T2 dx + - \ B(u+,u+)dt + - B(u+,u v )dt 



/I rTz pT 2 
<fi(glogg + l- g)\ Tl dx - - J B(u+,ip)dt + J Q(u+,u)dt. 



(iii) We will think of the LHS as the basic energy of this calculation, and the RHS 
as the terms still to be controlled. 

(6|3. Cutoff functions, control of the RHS and final goal 

In order to tackle the RHS and continue the proof of the lemma, we need to make 
a convenient choice of the sequence of cutoffs. Though only some simple bounds on 
the functions and their derivatives are used, a possible practical choice is as follows: 

(6.6) M%) = min{l + (|x| £ - 2) + , <p k (x)}, 7p k {x) = 7 - + ^ - U~ k , 
for some small e > and k = 1,2, . . . 
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Note that (fk > (fik-i- The following remark will be important: at points where 
ifk < 1 we have 



1 A-k 

Lp k = (fk-x + -4 . 



We also have 



inf Lp k = v?fc(0) > 1/2 for k > 1. 

Moreover, ipoo(x) < 1 precisely for \x\ < y/2 and ipi(x) < 1 for \x\ < 2. This 
means in particular that fk{x) = 1 + (\x\ £ — 2) + for \x\ > 2, k > 1. Moreover, 
V?oo(x) = 0| 2 + 14)/16 < 15/16 for |x| < 1. 

A more general version of the same construction takes 

(6.7) a = 1 __L + M!_i c -, 

with C possibly larger than 4. In that case 1 — ^Poo{x) > 1/4C for \x\ < 1. 

For the rest of this proof we write u\ = {u — <Pk) + > 0, u k = {u — ip k )~ < 0. Notice 
that the support of u\ is contained in the ball of radius 2 as a consequence of the 



assumption (5.2). 



We are ready to tackle the RHS of Identity (6.5) with this choice of test functions. 



One part will be controlled by a small multiple of the present energy, i.e., we will 



absorb it into the LHS of (6.5). The rest will be bounded above by a large multiple 
of \{u\ > 0}l (t ne notation |.| means the measure of the set). We recall our goal: if 
we do this, together with Sobolev inequality, we will get an iterative relation for the 
LHS energies 

A+i < C fc (A) 1+a 

that converges to zero as k — > oo, as desired, if the iteration is started from a small 
initial value Aq. 
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[6[4. Estimate of the remaining B term 

We start the process with £>2 = B(u^, ipk)- By inspecting the integral we easily get 

B 2 <jB(u+, u+) + -B*(<p k ,<p k ), 

for every 7 > 0, where B*(<p k , ifik) indicates that the integral is performed only on the 
set where either x or y belong to {u^ > 0}. That is, 



B * = J J ^ u t >0 ^ + x H>o}( y ^ K ( x ~ V)^k{x) - <p k (y)) . 

For 7 small, then 7-B(m^ ,u^) is absorbed into the LHS (into the energy). Now, using 
that 

\<Pk(x) - <Pk(y)\ < Cmin(l, \x - y\), 
and the size of the kernel K, we arrive to the estimate 

^<C|K>0}|<C4 2fc f {uUfdx 

J{u+>0} 

The last inequality follows by Chebyshev's inequality, since Uk-i > 4r k /2 whenever 
Ufa > 0. The obtained expression is good for our later purposes. 

[6[5. Analysis of the Q terms for < s < 1/2 



The last term in (6.5) also has a bilinear structure. Indeed, 



u t ( x ) — VL(x - y)[u(x) - u(y)]dxdy := Q(u^,u) 

= Q«A) + Q«^k) + Q«,%) = Qi + Q2 + Q 3 , 

but note that the "kernel" that is involved is not symmetric due to the presence of 
terms with (p k - The study of the contribution of each of the three terms is again split 
into the close-range and far-field interactions, represented by the integrals for x — y 
lying in a ball around the origin, or in its complement. In that sense we note that 
VL satisfies |VL| < c\x — y\ K(x,y). This will be used repeatedly. 
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Let us first tackle the integral Q\ in a ball of radius i] around the origin 



\x-y\<ij 



u+(x) 



c\x-y\K(x,y)\u£(x) -u+(y)\dxdy 
< 8c 2 / / (ul) 2 (x) \Vip k /(p k \ 2 \x -y\ 2 K(x,y)dxdy 

J J\x~y\<n 



'\x-y\<Ti 

K{x,y)\ul{x) -u+{y)\ 2 dxdy . 

'\x-y\<v 

The last integral is just the part of (l/A)B(u^,u^) integrated for \x — y\ < rj, so it 
can be absorbed by the LHS energy. The first integral is first integrated in y which 
is easy since 

K(x,y)\x-y\ 2 dy = 0(r] 2s 



\x-y\<?) 



Using this and also that |Vy?/V| < C, this first integral can be estimated then as 

< Crj 2s J \u + k f dx , 

an expression that can be left in the RHS or absorbed into the LHS if we take 77 small. 

• Let us compute the outer part of Q\ (for \x — y\ > rf). In this region, VL is 
integrable, so that 

|Q ou 'KW)| < C( V ) (| {utfdx + ( J v+dx)- 

which is also admissible, as we will see. Note that the last integral comes from the 
term in u^(x)u^(y). 

• Next, we treat the term 

0.2 = Q{ut,ipk) =11 4— - V)[fk(x) - ip k (y)]dxdy. 

JJVk 

Remember that u%(x) is compactly supported in a small ball. For \x — y\ < 4 we 
have \<pk(x) - <Pk(y)\ < C\x — so 

VL(x - y)[ipk(x) - Vk{y)\ 

is integrable in y and we are left with 

C [ ut(x)dx < C\{u + k > 0}| <C k [ (uUfdx, 

J J{u+>0} 
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which is also a good term. 

• For \x — y\ larger the calculation is more involved. We will use the fact that VL has 
mean value zero on spheres, Therefore, we write ifk(x)—ipk(y) = (ipk(x)—(pk(y—x)) + 
((fk(y — x) — <fk(y))- We observe that 

/ VL(x - y)((pk(x) - <p k (y - x)) dy= VL(z)(ip k (x) - (p k (z)) dz 

J J\z\>4 

is zero in the sense of principal value since (i) (<fk(x) — fkiz)) is a radial function of 
z, (ii) we have an antisymmetry property for VL; both facts imply the cancelation of 
the integral. The rest of the integral is 

J VL(x - y)((p k (y - x) - if k {y)) dy 

Since \VL(x — y)\ ~ \x — y\~( N+1 ~ 2s as \y\ — > oo and \ipk(y — x) — Pk(x) \ ~ C\x\ \y\ £ ~ l 
the whole integral is convergent if 2 > 2s + e, which is a smallness condition on e - 
Performing then the integral in x, we get the conclusion that 

\QW,v>k)\ < c|K > °}l < ck [ K-i) 2 dx, 

J{u+>0} 

as desired. 

• The last term to examine is 

(6.8) Q3 = QK,u k ) = -JJui(x)?^VL(x- y )u- k ( y )dxdy. 

For \x — y\ < rj small, we use that \VL(x — y)\ < C\x — y\ K(x — y) and then Q(u k , u k ) 
is bounded by a small fraction of B(u k ,u k ) (remember that this term had the right 
sign). We can therefore get this part absorbed by the LHS of the energy identity. 

• Finally, for \x — y\ > i], we have the worst convergence case. This is the only place 
where we use the restriction s < 1/2. We solve the difficulty of the integrability in y 
at infinity by taking e < e = 1 — 2s, so that integration first in y is bounded, since 
the term \VL(x — y)u k {y)\ is of the form 0(\y\ 2s+e ~ N ~ l ). We conclude that 

(6.9) |Q OMt «,%7)| < C|K > °>l < Ck I K + -i) 2 - 

J{u+>0} 
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Summary. Using (6.3), we obtain for < s < 1/2 and t± < t 2 < the following 
energy inequality: 



(6.10) 



/ 



t-2 



< 2 



<Pk 2 J tl 



B(u+,ut)dt 

-t 2 



Vk 



-dx + C 



2k 



ti J{u+>0} 



( u t-i) 2 dxdt . 



where C is a universal constant that only depends on s and the dimension, N. In the 
application to the iteration, the ti will be chosen in dependence of k. 



[6[6. Iteration and end of proof of Lemma 5.1 



This part is very similar to the one at the end of the boundedness proof in section 
We define the total energy function for the truncated solution as 



4.2 



(6.11) 



A k = sup l(u~l) 2 (t)dx+ / B(u~£,u~l)dt, 

T k <t<0j JT h 



where T k = —2(1 + 2~ k ), k = 0, 1, • • • . Notice that tp k lies between 1/2 and 1 at 
the points where is not zero. From (6.10) with k > 1, taking arbitrary values 
t 2 = t >T k and t\ = t' G pfe_i, T k ] we have 

(6.12) A k <A inf f(u+) 2 (t')dx + C 2k [ [ (n+_J 2 dxdt = I + 1 1 . 

i'epVxjy J J t , J{u+ > Q} 

Taking averages in t' we arrive at the inequality 



inf 



(u+) 2 (t')dx< 



< T 



T k -i 



T k - T k _ x } Tk _ x 



'ut)\t') dxdt'. 



{u+) 2 {t')dxdt' 



Observing that u^(x) > implies _ 1 (a^) > u^(x) +4 fc /2, we can realize that both, 
/ and 77, have the same flavour, and that in fact we have the estimate 



(6.13) 



A k <C k 



T fc _i-/{<_ 1 >4- fe /2} 



( u t-i) 2 dxdt . 
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for a possibly larger constant C. The next step it to modify the proof in section 4.2 



replacing the Sobolev inequality by the second part of Lemma 3.2 To that end, let 



p > 2 be the exponent corresponding to Sobolev's embedding theorem so that 

Take 9 = 2/p and define q = (1 — 6)2 + dp. Then 

«_i) 2 dx < 4( fc+1 ^- 2 ) / («+_!)« dx 



f 

< 4( fc+1 )(<?- 2 ) (J {u+^fdx) (J (u+^ydx 



Integration in time t along the interval [Tfc_i, 0] gives us from inequality (6.13) and 
the previous estimate a recurrence relation of the form 



Ak < C k (sup Tfc i < 4 < J (n+_ 1 ) 2 (t) dx^j ■ B(u+_v dt 
with t = 1 — 9 > and a possibly larger constant C . 

We need 5 to be very small to start the iteration so that the sequence Ak converges 
and then = 0, which means that u < 7]^ and this in turn implies that u < 7/8 
for |x| < 1. We thus get the result in the Lemma statement with p, = 1/8. 

Remark. A simple modification of (p,^ would give other values of p, £ (0,1/2), of 
course with a different estimate of the maximum allowed value for 5. The proof also 
shows that the time size T = 4 can be replaced by any other number and the iteration 
will work with a different value for S (and the same values for p and e). □ 



7 Modification of the energy calculation 

In the iterative process that we will consider below there will be situations in which the 
solution considered in a cylinder as above is bounded between two positive constants 
< Mi < Mi- We want to establish that a similar result holds and the relation 5-p 
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does not change much, which will essential in our iterations. The use of rescaling 
and the translation invariance in (x, t) allow to recover a solution defined in the 
standard domain T 4 which is the one chosen for all our calculations. But imposing 
the normalization < u < 1 asks for a vertical translation in u to adjust the lower 
level (on top of the usual scaling), and this leads to a modified equation with the 
following form 

(7.1) d t u = V • (D(u)VCu) 

where D has the form D{u) — d± + d 2 u. We will normalize so that d\ + d 2 = 1 (i. e., 
D(l) = 1) and we will have \d 2 \ < 1 (in practice, d 2 becomes small as the iterations 
advance). 

We re-do the energy calculations of the previous section but this time we use as test 
function rj = F(u V ip) — F(ip), where F is defined as 

f u 1 1 

J 1 D(s) d 2 



Note that rj > and rj = where u < (p. Then the LHS gives 

«2 J Tl J \ d 1 + d 2 tp 



hhs = i [ rjdtudxdt = -j- [ /log ( d l±^lll ) dtudxdt. 



Next, we write 



d 1 + d 2 (uy v ) =i+ fi= d 2 



di + d 2 ip Lp) di + d 2 ip ' 

and note that the function 

H(s) := ^-{s log(l + (3s) + - log(l + (3s) - s} 
d 2 (3 

is such that H'(s) = d 2 1 log(l + (3s), so that by integration in time we arrive at 

I lhs = J H{u+{x,T 2 ))dx- J H(u v p + (x,T 1 ))dx. 

We will need an estimate of this quantity. Since H(0) = H'(0) = and H"(s) = 
(3/(d 2 {l+/3s)), i.e., H"(u^) = l/(d 1 + d 2 (uV(p)), we have the estimate if(uj) ~ (mJ) 2 
at all points where u > ip, u is bouded above and ip is bounded below away from zero 
in the sense that 

(7.2) Cl «) 2 <//Kp + )<c 2 «) 2 . 
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and the constants go to 1/2 as d 2 — > (and d\ — > 1), since in the limit H"(s) = 1. 



On the other hand, on the RHS of the weak formulation, instead of Formula (6.4), 
we have the following: 



dt j (F(u V<p)- F(ip)) div [D(u(x)) VL(x - y) (u(x) - u(y))} dxdy 

-I a I { V "W v ^> 

./ Ju>tfi - 



-} -D(w(x)) V x L(x — y) (u(x) — u(y)) dxdy . 



ixp D ( u ( x )) D{ip{x)\ 
which we again split as I + II. In the present situation we take 

I = — J dt J J VmJ(x)V ' x L(x — y)[u(y) — u{x)]dxdy. 

After integrating by parts and symmetrizing, we get —I = (1/2) f B(u+u) dt, where 

<B(mJ,u). = J J (u+(x) - u+{y)) K(x - y){u{x) -u{y))dxdy. 

As before, we separate the expression into three integrals, using the splitting: u = 
■u+ + (p + u~. The rest of the integral in the RHS takes the form 

11= I dt H Z ^^ Vip(x)VL(x - y)(u(x) - u(y)) dxdy. 

J JJu(x)><p(x) D \ l P) 

• We note that when D(u) := d\ + d2U with d\ + c?2 = 1, we have F(s) = 
(l/d 2 ) log(di + d 2 u) and II = / Q{ut ) u) dt, with 



Q{u+u)=d 2 u+{x) 



Vip(x) 



-VL(x — y)(u(x) — u(y)) dxdy . 



^ J D(<p(x)) 

which looks like Q of previous section but for an interesting small factor, d 2 - 
Finally, the energy inequality takes the form 



(7.3) 



#«) dx 



r 2 

Ti 



+ C r B(u+,u+)dt + C f 2 B{u%,u~)dt 



<C 2 / B(u+<p)dt+ / Q(u+u)dt 



T-2 



to be compared with identity (6.5). 

Repeating the rest of the steps of the previous section offers no novelties and we 
arrive at a similar result 
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Lemma 7.1 The result of Lemma |5 . 1 1 holds for the weak energy solutions of Equation 
(7.1 ) under the assumptions D{u) — d\ + d 2 u, di, d% > 0, di + d 2 = 1, d\ > 1/2, and 
he result holds with same constants fi, 5 for different values of d\,d 2 . 



8 Pulling up from zero. Proof of Lemma 5.2 



Here, we are trying to pull the equation uniformly away from zero under some as- 
sumption on the size of a level set. The precise assumptions are 

< u < max{l, \x\ e - 1} in M w x[-4,0], 

and \{u < 1/2} fl T 4 | < 5|r 4 |. The desired conclusion is then that u\ T2 > /i > 0, for 
an appropriate 5 > to be chosen. 

The tools are energy inequalities and integration by parts. In order to deduce the 
necessary energy inequalities we recall the bilinear forms and integration by parts 
formulas of Section [6] We use the same notation for the kernels K and L. 

[8[l. Local energy inequality. The basic calculation is as follows. We take a posi- 
tive smooth function ip and use rj = log(v) as a test function in the weak formulation, 
where v (x, t) is defined as v = u/ip if u < if and v = 1 otherwise. In other words, 

(8.1) v = ^M, or v-l= { ^^:=^. 

if if cp 

Note that < v < 1. We will use the notations: w+ = {u — f) + , u~ = {u — f)~ , so 
that one has u = f + u+ + u~ . Next, we put 

h(v)=[vlog(v)-(v-l)], Vw>0, 

This function takes as minimum h(l) = and h(v) > for all v > 0, v ^ 1 with 
h(0) = 1. Moreover, h"(v) = 1/v so that for < v < 1 we have h"(v) > 1 and 

1 (u~) 2 

(8.2) h(v) = [v \og(v) + 1 — f ] > -(v — l) 2 = . 
We also have the inequality 

(8.3) h(v) <(v- l) 2 = K -^-. 

f^ 
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This is used below. Next, we have the following calculation for the choice of function 
v = v(x,t) made in ( |8.1 ) 

^ iph(v)dx= [\og(v)ipv t = [ \og{v)(pv t = ( \og{v)u t 



dt 



u<ip 



since \og{v) = for u > tp. Using the weak formulation of the equation, we now get 



U) 



iph(v(t 2 )) dx + 



<*p h(v(ti)) dx. 



j V[log(t>)] uVCudxdt = j 
Let us work out the meaning of the middle term (energy term). First, we have 

/V[log(v)l - uVCudx = [ u— — - VCudx = [ (pV(u~/(p) - VCudx 



5.5) 



Viiy ■ VCu dx — 



{u<<p} 

^V<^ • VCudx. 

{u< V } V 



Clearly, using u = ip + u+ + u v = (u V cp) + u v , we get 

/ Vw~ • VCudx = B(u~,u) = B{u~,u~)+B{u~,u+)+B(u~,(p). 

The first term is the one we want to keep to complete the expression of the energy. The 
second term B(u~,u+) is also positive in view of formula (3.1) since both functions 
have opposite signs and disjoint supports, cf. Corollary [3TT Hence, this term could 
be discarded, but it will turn out to be useful as we have already seen. 



The remaining term B(u~, <p) and the last integral in (8.5) are delicate and it is there 
that we have to make an argument with a careful choice of test function. Summing 
up, we have 



V [ \og(v)} ■ uVCu dx > B(u-,u~) + B{u~, <p) + B{u~ , u; 
Putting 



u<ip 



<P 



-Vy? ■ V Cu dx . 



Q(w,u) = JJ 



W[X) 



Vp{x) 
<p(x) 



VL(x — y)[u(y) — u(x)]dxdy, 



we arrive from the above and from (8.2) and ( |8.3 ) to the basic energy inequality: 
1 



(8.6) 



< 



o ^(u-(t 2 )Ydx + 

- (u~{ti)) 2 dx 

ip 



rt2 rt2 

/ B{u~,u~)dt+ / B(u~,u^)dt 
Jti Jti 



t-2 



tl 



t-2 



B{u-<p)dt+ / Q(u~,u)dt 



ti 
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The local energy, in the time interval [Ti,T2], is now defined as 



(8.7) 



E v (u) := sup 

Ti<t<T 2 




— [u ] 2 dx + 



B(u ,u )dt. 



Note that the test function log v is negative but we are interested in u being a super- 
solution, and this is the case for instance if we truncate it by 2: u = u A 2. So, we do 
not need to worry about growth at infinity. 

|8}2. The iterative process 

(i) At this point, following De Giorgi's idea we want to obtain an iterative relation 
playing the energy inequality against the Sobolev embedding using a convenient se- 
quence of cutoff functions. We consider a series of cutoffs with dyadic separation, so 
that 



holds in the support of (pk+i for some fixed < a < 1. Moreover, the ipk converge 
to /io from above in -62(0). We will use as a test function r\ = logv with v = u/<p>k 
if u < (fk and v — 1 otherwise. We arrive at the energy inequality at the end of last 
subsections with ip k instead of (p. 

(ii) We propose a concrete construction of the cutoffs (pk- All functions <fk{ x ) are 
symmetric, nonnegative, non-increasing and compactly supported and the sequence 
is nonincreasing with k. The basic profile (po is a kind of rounded mesa: 



for some m > 2 to be chosen later. 

- To construct <pk we first rescale the graph of tpo from the interval [3, 4] = [2 + 1,2 + 2] 
to the interval [2 + 2~ fe ,2 + 2 _fe+1 ], and extend by fi + (1 - /io)2~ fc inside the ball 
B 2+ 2-k. Then, ipk has the following properties: 



-k 




= \[A - \x\] m when 7/2 < \x\ < 4, 
= 1 for \x\ < 3, 

a C 2 decreasing radial function for 3 < |x| < 4, 





(d) 




|v^|M<cV, /m , 

<Pk-i -¥»*>(!- Ato)2~ fc , 
¥k -> VoXb 2 as k -)■ 00, 



in the support of <^fc, 



so that a = 1 — fiQ. 
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1 




(iii) Sobolev embedding. The embedding that we need is a variant of the one used 
in preceding sections. We have seen that B(u~, u~) is the square H r norm of u~, and 
in this way it controls an L p norm of u~ for a given p > 1 that we have calculated. 
Coupling this with the energy term 



sup [ —[u ] 2 dx . 

<t<T 2 J 



Ti 

we can control for some q > 2 and 6 > 

{£ J<p^\u-\*dxdt) 2lq . 

(Since the cutoff <p goes to zero, the term ip~ e will be important in controlling the 
term containing Vifi/f). The proof is as follows: Write q as a convex combination, 
q = 26+p(l-6). Then, 



J tp- e \u-\« dx < ( J V - l \u-\ 2 dx) 6 { J \u-\vdx) 



(1-9) 



We make the choice 6 = (p — 2)/p with p as the q in Lemma 3.2 Then, 1 + 6 = q/2 
and 

,T 2 , N 2/9 



y ip- e \u~\ q dxdtj < e v . 



From now and on we will write and Ek to denote and respectively. 
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[8|3. Iterated energies. In the left of (8.6) we truncate in time along an increasing 
sequence Tk — > —2 and get, for every Tk-i < t\ < Tk < ti < 0, 



(8.8) 



2ip h 



(u k {t 2 )fdx 



t2 



*2 



B(u k ,u k )dt+ I B(u h ,uJ)dt 



< / — (u k (ti)) dx— I B(u k ,ip k )dt - / Q(u k ,u)dt. 



t-2 



|8}4. Analysis of the RHS 

We now examine the terms left in the RHS of the energy inequality (8.8) that we 
will call I, II, and III for convenience of reference. Our purpose is to show that 
these extra terms are either bounded by a small multiple of the LHS, so that it will 
be absorbed by the LHS, or they are controlled by a term of the form 

C k J 2 J ip- e \u-\ 2 dxdt. 

Estimating I. As initial times Tk we will choose the time T k = — (2 + 2~ k ) and 
t\ G [Tfc_i,Tfc] is the point where the value of 

inf / — (u k (t)) 2 dx 

T fe _i<i<T A J if k 

is attained. In this way we have 



-2-2 J - fc J ¥k 

Since \{u — <fk)~/fk\ < 1 we have [(u — <£k)~} 2 /<£k < \{u — <Pk)~\-> an d we can bound 
/by 



C k i 

Estimate of II. About this term we have 



X{u~<0}- 



| B(u k , ip k ) | < ^ [B(u k ,u k ) + B* (ip k , <p k 
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where the star indicates that we can restrict the domain of integration defining 
B((p k ,ip k ) to the points where u k (x) 7^ or u k (y) ^ 0. So this last term can be 
replaced by the better expression 

^// - ( Pk(y)] 2 K(x,y)x(x,y)dxdy, 

where x is the characteristic function of the set of points (x, y) where either x G 
supp u k or y G supp u k . From the Lipschitz character of ip k and that \V(p k \ < C2 fe 
we have for fixed x 

J [ip k (x) - <p k (y)] 2 K(x, y) dy < C k . 
Therefore, the above term may be bounded by 

C k J X {u - <0} dx. 

Hence, 

n <l[ 2 ®^~ k ,u- k )dt + C k [ 2 J X K -<o } dxdt. 
The first term is absorbed by the local energy in the LHS. 

Estimate of III. The remaining part of the proof is devoted to the estimate of the 
term 

VMX) VL(x-y)(u(y)-u(x))dxdy 



Write, by the usual splitting of u, 

Q{u- k ,u) = Q{u k ,u k ) + Q(u k ,u+) + Q(u k , Vk ) = Q 1 + Q 2 + Q 3 . 

<2i : Noting that 

\VL(x,y)\ < \x-y\K(x 1 y) = 0(\x-y\- { - N+l -^) 1 



we have from Holder's 

\Qi\ = |Q(«fc,Wfc) 



jju k {x) 



VL(x - y)[u k (y) - u k (x)]dxdy 



< 



< 



u k [X] 



Vip k (x) 



1 



\x-y 



N-2s 



dxdy 



1/2 
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<4CV J V - k 2lm {u- k Ydx+\B(u- k ,u- k ), 



where we have used here the properties of the functions f k . Integration in time gives 
that the second term is absorbed by the local energy expression in the LHS and the 
first is admissible for our purposes. 

Final step. To study <2 2 and Q 3 we consider a smooth decomposition of the kernel 

VL = VL^, + VLi_ f 
where 1 — ip is a bump function supported around the origin. We get 

Q 2 + Q 3 = J J u k (x)^^VL 1 ^(x-y)[ut(y) + My)]dxdy+ 

+ 11 ^) V () ^L^{x - y)[u(y)]dxdy =: Q 2 + Q 3 . 

In other words, in the compact part of the support of 1 — ip we keep the expression of 
u as u k + ip k + u k , while outside, where u k = and ip k grows, we just keep the term 



/ 



VL^(x - y)u{y)dy. 



• This term has to be controlled in a different way (through the change of coordinates 
if s > 1/2). As a consequence, the integrations in y for Q 2 and Q 3 are convergent 
and we are left with estimating 

/ ^ ^^(x) j Vk 2 ' m ( U ~k) 2dx + jx {u -<0} dx - 

The integration in time produces again two admissible terms. 

• Finally, for the term 

u k (x ^^f^ VL^x - y)[ul(y)]dxdy, 



we use the good term 



B ( u k . u t ) = J J u k (x)K(x - y) [u+(y)]dxdy, 
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left in the energy inequality, to absorb the integral whenever 

Vy? fc (x) 



<fk{x) 



\x — y\< r\. 



In the complement, that is when 



\x — y\> rj. 



we use that </? fc 1//m > |Vy2 fc |/<y2 fc and integrate in y: 



VLi^{x-y)[ul{y)]dy 



< 



4 n-1 



-dr ~ max ((J, rjip^ s ^ m ^j _ 



In any case, we end up in the worst of the cases with an expression of the form 



—s/m I — I ; 



that we control as before. 



Inserting estimates II and III in (|8.8|) we get that 

/ c 

(8.9) 



1 



< / — {u k {t 1 )) 2 dx + C k 



.) „ k K (*2)) da; + ^ £(u fc , w fe ) dt 

/ ( Pk(u^) 2 dxdt+ j 
ti J Jti 



X{ u - <o}d%dt 



whenever T^_i < ^ < < t 2 < 0. In particular, 
(8.10) 



su PT fc <t 2 <o / ' ^-{u k {t 2 )) 2 dx+ ! B{u k ,u k )dt< inf [ ^- {u k {t 1 )) 2 dx 



+c k 



< c k 



/ ¥k~ 2/m (Uk ) 2 dxdt+ / 



Or r0 

tp k 2/m (u k ) 2 dxdt + 

T k -i J JT k _ 1 



X {u -<o} dxdt 



X{ u - <0 }dxdt 
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where the last inequality comes from estimate J. In fact, the same argument given for 
this estimate shows that if m > 2 then the leading term in the last expression is the 
second one. Thus, if we define 

(8.11) A k ■= sup (u k ) 2 dx + j B(u k ,u k )dt, 

T k <t<0j JT k 

then, (8.10) gives 



o 



(8.12) A k <C k / X{u - <0} dxdt. 

Using that ipk-i > fk + a2~ fc in the support of <f k we get than 1 > \u k \ + a2~ 
in {u k < 0}. Therefore, 



(8.13) A k < C {1+q)k [ [ \u k ^\ q dxdt. 



The rest follows the same argument as in Lemma 5T let p > 2 be the exponent 
corresponding to Sobolev's embedding theorem so that 



J lu^dxKCiBiu^u^f 2 . 
Take 9 = (p — 2)/p and define q = 62 + (1 — 6)p. Then, 

A k <C^ k £ (/l^-xl 2 ^ B{u- k _ v ul_ x )dt 
^ 8 ' U " > < & ^sup Tfc _ 1 < t < J (Wfe_!) 2 (t) dxj jf SCWfc-i^fc-i) dt 



< C k A{_ x A k ^ = C k Alt{, 



This completes the proof of Lemma [5. 2| for the critical case since, by hypothesis, 

A = y° J [(m ~^ o) ]2 dxd< < jj^ xiu-^o} < \{u < 1} n r 4 | < <5|r 4 |, 

and S can be chosen as small as we need. 

The noncritical case of Lemma 15.21 This is the case where the solution lies 
between, say, M and M + 1 for some M > 0. The proof is similar except that for the 
lower estimate we may already assume that w|r 4 > p> > and then all the involved 
cutoffs satisfy ip k > \i so that V(p k /(p k is a smooth bounded function. 
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9 The lemma on intermediate values 



We still need a main ingredient before we attack the regularity issue by means of a 
suitable iteration. Indeed, we have to improve Lemma 5T by showing that, in order 
to get a uniform reduction of the maximum in a smaller ball it is not necessary to ask 
that u < 1/2 "most of the time" in T 4 , but only "some of the time". This is precisely 



stated in Lemma 5.3 



The proof of this result uses De Giorgi's idea of loss of mass at intermediate levels, 
which he applied in an elliptic context. In the present parabolic setting we will follow 
the idea of the proof of the similar result that was carried out in the papers [10], [6] 
in the linear framework. We give next the detailed statement of the "intermediate 
lemma" and its proof. We start by selecting a cutoff function of the form <p = 
1 + + where we choose 

F{x) = sup{-l,inf{0, \x\ 2 - 9}} 

Note that F < is Lipschitz, compactly supported in -63(0) and F — — 1 in B 2 (0). 
Moreover, for < A < 1/3 we define 

if>x(x) = ((\x\ - \- 1/u Y - 1)+ for |x| > \- 1/u , 

and ip\(x) = otherwise. This represents a "wing" that starts far away when A is 
small, as will be the case. A convenient value of v G (0, 1) is needed and it will be 
determined later. We also define 

<pi = 1 + + XF, V2 = 1 + ^j x + A 2 F, 

and put (fo = <p, so that ipo < ipi < ip 2 < 1 in the ball -64(0). 

Lemma 9.1 There exist small constants p, 7 > depending only on N and s, and 
A G (0, 1), depending only on N, s, p, 6 such that for any solution u defined in 
with 

u(x, t) < 1 + ip\ in S4, 

with X < Ao, and also such that \{u < (fo} R ((—4, —2) x B\)\ > p, then we have the 
following implication: if 

\{u> V2 }n((-2,0)xR N )\>5, 

then 

(9.1) |{v? < u < <p 2 } n ((-3, 0) x R N )\ > 7. 



38 



The last line asserts that under the stated assumptions the measure of the interme- 
diate level cannot be small 

Proof, (i) In our context we start from the energy estimates we have obtained during 



the proof of Lemma 5T We have to arrive at a "correct form" of the energy inequality. 
We recover the energy calculation (6.5) with <p equal to the "intermediate cutoff" ipi, 
i. e., 



(9.2) 



(fiH(g) dx 



T 2 
Ti 



+ 



T 2 



B((u - <pi)+, (u - <p 1 ) + )dt+ 



1 r 2 1 f 12 



Tz 



Qdt, 



Ti 



with g = l + (u — (pi) + /ipi and H(g) = g logg-(g-l). We want to get now estimates 
on the RHS that show that all the terms are either absorbed by the LHS or can be 
estimated above by C\ 2 , where C is a fixed constant. With this we will arrive finally 
to the expression 



r t 2 r T 2 

/ <Pi(g \ogg + l- g)dx +c B{(u - <fi)+, (u - <pi)+) dt 
+c / 2 B{{u - ^i)+, (u - ^)_) dt < CA 2 (T 2 - Ti) . 



(9.3) 



This would complete the preparatory step of the lemma. 



ii) The verification that (9.3) holds is as follows. For simplicity we write (u — ip 



u{ and {u — (pi)- = u 1 . Repeating a bit some arguments, to prove (9.3) we observe 
that 

|£«y?i)| < jB{u+,ul) + -B^cpt^i) , 



7 



for every 7 > 0, where 



Xb 3 {x)K(x - y){ifi{x) - ipi{y)fdydx. 



For 7 small, then ^B^^Ui) is absorbed into the LHS of (9.2). For the second term 
we have 

£*Wi) <4(A 2 £*(F,F) + £*(^ A )) 
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Since the function F is Lipschitz, compactly supported and \x — y\ 2 K(x — y\ is 
locally integrable, we have B*(F, F) < C. Now, for A small one has X~ x l v > 4 so that 
ip\(x) = if \x\ < 3. Therefore, using that ip\(y) < \y\ u always, we get 

S*(Va, V>a) = [ [ K{x- y){Mv)) 2 dydx 

JB 3 J\y\>\-V" 

JB-i J\y\>\-^l» \X — y\ Jb 3 J\y\>\-i-/v \y\ 

It suffices to take 2 ~ 2s ~ 2v > 2, that is, v < 2 -^ 1 i in order to conclude that 

B*( Vl , Vl )<0(\ 2 ). 
We now consider the terms involving Q: 

Q(ut,u) = Q{u+,uf) + Q{ut,vi) + s«,«r)- 

The main ingredients are the following estimates: 

u t < -^1-^1 > since < u < 1 + if)\. In particular, 
wf < Axs 3 

|VL(x-y)| ~ |x — 3/|iC(ic — j/) ~ [^py+r=2i , integrable at infinity for s < 1/2. 
\V(pi(x)\/ipi(x) < C, Vx, with C independent of A, and 
\Vipi(x)\/(pi(x) < CX, in the support of uf. 
Hence, 



|Q«,OI = 

<Cr/i3«, O + CA /Y [uf(x) 2 + uf(x)uf(y)]\VL(x-y)\dxdy 

J J \x-y\>rj 

< Ci]B(ut, uf) + CX [ ut(x) 2 dx + CX( f ut(x) dx 



The first term in the last line of the formula goes to the LHS of (9.2) for rj small and 
the other two are just (9 (A 3 ). Also, 

Q(uf,< fl ) = XQ(ut,F) + Q(ut, V>a) 
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Since 

\Q(uf,F)\<CX J uf(x) J \VL(x -y)(F(x) - F(y))\dydx <CX J u+(x)dx, 
the first term above is easily seen to be of order 0(X 3 ). For the second we have 
\Q(utM\ <CX [ uf(x) [ \VL(x-y)(Mx)-My))\ dydx 

dy dx, 



< CX / ut(x) 



i„i>a-v. \y\ N+1 ~ 



2s 



where we have used that ip\(x) = in the suport of and ip\{y) < \y\ u ■ From this 
we get 

\Q(utM\ < CX 1+1= ^ J ut(x)dx. 
Observing that v < 1 — 2s from the previous condition on v , we conclude that 

Q(u+^ x )<0(X 2 ). 

Finally, 



ie(«r,«r)i = 



tpi(x) 



-VL(x — y){u l (x) — u x (y)) dxdy 



u{{x) 



\x-y\>r) 



V^i(x) 



ip x [x) 



\VL(x - y)u t (y)\dxdy. 



The term CriB{u{ ,u 1 ) is absorbed into the LHS of (9.2) for rj small. For the other 
we use the fact that \ui \ < (fx — 1 + XF + Now, 



u t( x ) 

"|V|>A-V" 

< CX [ uf{x) 



ipx(x) 



\VL(x-y)(l + i> x (y))\dxdy 



i + M' 



|y|>A-V 



- \y\ 



N+l-2s 



dydx<CX 1+ v j u^(x)dx 



whereas, 



u t( x ) 



<Cxjuf{x)J ^ j- - ^ N+1 _ 2s dydx < CX J uf(x) dx. 



Vy>i(z) 



(px{x) 



\VL{x-y){l + F{y))\dxdy 
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In both cases we get 0(X 2 ) for an appropriate v, which proves (9.3) as wanted. 



(iii) Estimate ( 9.3[ ) shows a property of one-sided Lipschitz continuity in time from 
above for the space integral of the LHS. On the other hand, we also have the inequality 



(p 1 H(g)(t)dx < / ip x (g - l) 2 (t)dx < / uf(x) 2 dx < CX 



Vt. 



Inserting this into the energy inequality (9.3) we arrive at 

/ 2 B(uf, ul)dt < CX 2 , for - 4 < T x < T 2 < 0. 
We are now in the position to apply Steps 2 and 3 of the proof of Section 4 of [6] with 



only technical changes that we will omit, and we will thus get Lemma 9.1 Note that 
the quadratic (or at least superlinear) dependence on A in Formula (9.3) is absolutely 
necessary for the proof to work. □ 



10 Oscillation reduction result 



We are now ready to prove the strong oscillation reduction result, Lemma [573] , Though 
the general idea of the argument follows closely Section 5 of [6] , an interesting modifi- 
cation is needed to accomodate the nonlinearity of the equation: in doing the scaling 
of the solutions in each iteration we will now find solutions of a family of related 
equations, and we have to use the modified estimates of Section [7], that hold for that 
family. At some step we will have to apply the results of the previous section to such 
a family. This is an easy further verification that we will leave to the reader. 



Instead of Lemma |5J3 we will prove the following version that will be the one used 



later in the regularity argument. Let us fix some notation. For A as in the previous 
section, we define for any e > 

^,A(x) = ((|x|-l/A 4 / s r-l)+ if M>A" 4 / S , 
and zero otherwise. 

Lemma 10.1 Given p > there exist e > and /ii such that for any solution of the 
FPME in (-4, 0) x R N satisfying 

(10.1) 0<M<1+^ £ ,A, 
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and assuming that 

(10.2) \{u<<p }n(-4,-2)xB 1 \> p, 
then we have 

(10.3) sup u < 1 — fix. 

(-1,0) xBi 

Note: In the next section we will take p equal to 5o as denned after the statement of 
Lemma 5.2. 



Proof, (i) Consider jo — | ( — 3, 0) x -B3I/7, with the value of 7 given in Lemma 9.1 
We fix e > small enough so that 

^±<(M-/*-i) + 

for all x. We may take e = (s/4)A 2 - 70 for instance. For j < jo we consider the sequence 
of functions defined iteratively by 

(10.4) u j+1 (x,t) = ^(u 3 (x,t) - (1 - A 2 )) , 

starting from uo{x,t) = u(x,t), the solution of the FPME under consideration. By 
induction we assume that 

{Uj)+{x,t) < 1 + ^e,x{x) , (x,t) < (-3,0) x R N . 

So for j < jo we have Uj < 1 + We recall that A is fixed and small. 

(ii) Next, we have to check the equation satisfied by Uj, that turns out to be 

(10.5) u t = V(L' i (n)V£(u)) 

with a diffusion coefficient Dj defined inductively by the rule 

(10.6) D j+1 (s) =D J (\ 2 s + (l-\ 2 )), 

so that Dj{uj) = Dj + i{uj + i) holds. These are the type of diffusion coefficients for 
which we have proved the modified estimates of Section [7| We have to make the 
observation that, with the notations of that section, we have g^Wj+i = X 2 uj + i = 
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Uj — (1 — A 2 ) so that the integrals in Q do not change their form, in particular the 



formulas of the type (6.8), that depend on the integral 



VL(x-y)(u j+1 (y)-<p k (y)) dy 



are estimated by the same constants. Note that J RN VI(x — y) C dy = by antisym- 
metry of VL. 

(hi) By construction, the measure \{uj < ip } fl (—4,0) x Bi\ is increasing; hence, it 



is bigger than p for every j. We can apply Lemma |9TT inductively to Uj. As long as 
\{u k > if 2 } n ((-2,0) x R N )\ > 5 for k = 1,2, ••• ,j + 1, we have 

|{<^o < u j+1 < ^ 2 } n (-3, 0) xl w }|> 7. 

Therefore, 

\{u j+1 > V2 }\ < \{u j+1 > f }\ -7 < |{«i > <^ 2 }| - 7 < I (-3,0) x B 3 \-jj. 
This cannot be true up to jo, hence there must exist j < jo such that 

IK > ^2} n ((-2,0) x r n )\ < 5. 

We can then apply the first lemma to the next step, Uj + \. Indeed, 
u j+ i < 1 + ipx < 1 + ipi on (-3, 0) x 

and 

|K+i > l/2}fl((-2,0)xB 2 )| < |K+i>v?o}n((-2,0) x B 2 )\ 
< |K >^ 2 }n((-2,0) xR N )\<5, 



and then Lemma 5.1| implies that 

Uj+i < 1 — /i on (-1,0) x Bi. 
This gives the result with [i\ = A 2 - 70 //. □ 
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11 End of proof of regularity for s < 1/2 



The whole technical machinery is in place, and we are ready to prove Holder regularity 
by means of the iterative process outlined in Section [5| We take any point P = 
(x ,t ) G M. N x (0,oo) and prove that u is C a around P with an exponent that 
depends only on the parameters N and s of the equation, and a Holder constant that 
depends also on the L°° norm of u and a lower bound of to. 

We start with some reductions. There is no loss of generality in assuming that u is 
bounded in the cylinder Q since we know by Theorem 4.1| that u is bounded in any 



strip of the form Mr x (r, oo) with r > 0. Moreover, by scaling we may assume that 
to > 4. It will be then convenient to make a space-time translation and put Pq = (0, 0) 
assuming that the domain of definition of u contains the strip S4 = M. N x [—4, 0]. By 
scaling we may assume also that < u(x, t) < 1 in 5*4. 

Consider now a positive constant K < 1/4 such that the growth of the wings is 
controlled as follows: 

(H-1) l-{^/2) ^ £{Kx) ~ ^ A ' e(x) - 

The coefficient K depends only on A, fii and e > 0. The parameters are as in the last 
section. The iteration that we will perform offers two possibilities. 

• Alternative 1. Regularity at a degenerate point. Suppose that we can apply Lemma 



10.1 repeatedly because the lowering of the oscillation may be assumed to happen 
always from above. We consider then the sequence of functions defined in the strip 
S 3 = R N x (-4, 0) by 

(11.2) Ui+l(x , t) = —l— Uj(Kx , Klt) , K-J^y 

Note that this time the u/s are all of them solutions of the same equation. According 



to the running assumption, and using (11.1), we can apply Lemma 10.1 at every step 



so that we have Uj(x } t) < 1 — //1 in the cylinder Qi = B\ x (— 1, 0) for every j > 1. 
In view of the scaling ( |11.2 ), this implies Holder regularity around the point (0,0), 



where the solution necessarily takes the degenerate value u = in a continuous way. 
• Alternative 2. Regularity at points of positivity. After some steps of the iteration 



the assumption on the measure of the set {uj > 1/2} made in Lemma 10.1 fails. 



Then, we are in the situation where the oscillation is reduced from below thanks to 
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Lemma 5^2, which pulls the solution uniformly up from zero in a smaller cylinder. 
Then the equation is no longer degenerate, because after that step we have 



< fjf < Uj(x, t) < 1, 



in the cylinder B\ x (—1,0). Scaling the situation we will be in the conditions of 
the equation with diffusivity D{u) mentioned above, to which we apply either the 
modification of Lemma 110.11 or the modification of Lemma 15.21 In fact, 



we can 



apply the modification of Lemma 10.1 both from above and from below since the 
degeneracy has disappeared. In this way, we obtain Holder regularity at a point Pq 
where u(Pq) > 0. 



12 C a regularity for s > 1/2 



Here we want to prove the regularity result of Theorem 5.1 for 1/2 < s < 2. If we try 



to re-do the proof of the technical oscillation lemmas listed in Section [5j we find a 
convergence problem that is already apparent in Lemma |5.1| Indeed, the bulk of the 



proof contained in Section [6] works without modification, and a problem was found 
only in the last estimate of Subsection [6}4, regarding integral Q3 in an outer region, 
where the decay of the factors that we have been examining does not provide for a 
uniform convergence of the integrals if s > 1/2. 



12.1 Analysis of the difficulty 

A possible solution is to make use of the known fact that u(x, t) is an L 1 function in 
x uniformly in t, in order to bound the integral of the y terms in Q 3 , with integrand 
VL(x — y)u^(x), since L{x — y) is bounded for large \y\. Indeed, we can put together 
the different parts of Q in the outer domain, and after dropping the contribution of 
the term in u(x) since it is zero by the antisymmetry property of VL, we get 

u t( x ) — ^-VL(x — y) u{y) dxdy 
< C|K-,i)||i ( [ u+(x)dx) < C k \\u(;t)\\i ( / «-i(x)?dx 



This solves the problem in one application of Lemma |5.1[ but then 5 and /i would 
depend on llu^lK. However, to obtain the C a regularity result we have seen in the 
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preceding Section 11 that we need to iterate this and the other oscillation lemmas, 
we want to rescale and repeat, and then a difficulty re-appears, because we will keep 
expanding the u and the x and therefore expanding the integral at every step, so the 
constants will be ruined in the iteration. We need a way to control such behaviour. 

It will be convenient to examine the whole part of Q(u + ,u) that contains the diffi- 
culty, i.e., 

(12.1) [[ U + {x) YvM VL{Xj y )u{y)dxdy . 

JJ Vk(x) 

As we dilate and repeat in the iteration scheme, the term VL(x — y)u(y) also starts 
to build up as y tends to infinity. On the other hand, the integrability in y at infinity 
is lost if s > 1/2, since in that case VL decays like 

|VL| ~ \y\-( N -^ 

and this is not good enough. However, the good news is that 



V • (VL) ~ \y 



-(JV+2-2s) 



(valid for all second derivatives) which is integrable as \y\ — > oo. Noting that VL is 
integrable for y ~ if s > 1/2, we conclude that 

V(x,t) := (VL(x,y,t))* y u(y,t) 

has bounded Holder seminorm. Therefore, it would be enough to control it at just a 
point, for instance at x = (for an interval of times). 



12.2 The transport approach for the case s > 1/2 

The technical way to make use of the last observation is to perform by a change of 
coordinates x' = x — S(t) that introduces a transport term to counter the difficult 
term J VL{y) ^(y) u{y) dy. To be precise, we define 

(12.2) S(t) = J v(t) dt, v(t) = J VL (y) u(y, t) dy , 

and we observe that \v{t)\ depends also on u, and that \v(t)\ < C < oo since u(y,t) 
is in Ly uniformly in t. Indeed, the value of v(t) = v(t; u) is only controlled by the 
integral of u is space (i. e., the mass of u(t)). Next, we introduce change of variables 

(x,t) h-> (x\f) := (x - S(t),t) , 
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and we write the equation for u with respect to the new variables, u(x,t) = u(x',t'). 
The RHS does not change since we are performing a space translation for fixed time. 
However, the time derivative in the LHS transforms as follows: 

u t (x,t) = u t '{x' + S(t'),t') + (Vu)S'(t') = u tl +v-V x >u. 

The new term is what we are aiming at. The equation takes the convection-diffusion 
form 

(12.3) u t + v- V x ,u = V(uVCu). 

In the sequel we will write t for t' and u instead of u without fear of confusion. The 
space variable is still written x' . Next, we pass the term v ■ Vu to the RHS and 
multiply by log((w/y?) V 1), as we did in Section|6j to obtain the energy formula. We 
observe that in this case the RHS contains an extra term of the form 

V \og((u/<f) V 1) v(t)u(x', t) dx'dt . 

This integral must be computed only in the region where u > ip, and in that case 
(u/(p) V 1 = u/ip = 1 + («fc/Vfc)> so that 



1= ! dt I Vulv{t)dx'- I dt I v(t)dx' = h - h 

The first integral vanishes, and the second is precisely the troublesome term: 
dt / dx(ul(x')/(p k (x'))V(p h {x') I VL(y')u(y',t)dy' , 

Ju><fi k J 

that disappears in this way from the calculation. Alternatively, the disappearance of 
the bad term in the energy calculation in the new variables can be easily seen if we 
write the equation for u(x',t) in the more symmetrical form 

(12.4) u t = S7 X> U(j {VL(y - x') - VL(y)}(u(y) - u(x')) dyVj , 

to be interpreted in the same weak form, or weak energy form, that we used for 



u(x,t). In any case, this allows to prove Lemma 5.1 also for s G (1/2, 1) if we work 
in the new coordinates, and the constants involved in the result do not depend on 
the L 1 norm of the solution. The price to pay is that the slope of the distorted space 
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variables does depend on the w-mtegral. So, in the first step of the iteration process 
we have shown how to transfer the difficulty from a numerical term to a geometrical 
distortion. 

In order to sum up the result, let us introduce the bound M = 1 V swp t>0 v(t), that 
depends only on u via the norm sup t \\u(-, t)\\i. 



Lemma 12.1 Let 1/2 < s < 1 and let u be a solution of the FPME under the 
assumptions of Lemma 5.1, Let us perform the above change of variables so that 
is defined in smaller cylinder Ql where L = 4/(M + 1). Then the result of 



u[x 



Lemma 5.1 is true for u, with conclusion holding in a smaller cylinder Qx/m, & ma U 
depend also on M. 



Thanks to 12.4, it is then immediate to see that the modified Lemma 7.1, as well 



as the pull-up Lemma 5.3 are also true if stated in the form that we have used for 



Lemma 12.1 A bit more of attention to the details will show that the stronger 



reduction Lemma 10.1 also holds, since the iterations do not change the scaling in 
space and time. 



12.3 Analysis of the transport term in the final iteration 



When we try to perform again the iteration procedure of Section 11, one of the 
alternatives is repeated scaling around a degenerate point. In that case the iterations 
take the form 

(12.5) ^ +1 ( M ) = _L_ u ,(Kx, Kl t), K X = i 5 (V4) ■ 

that we may sum up as 

Ui(xi, tx) = A u(x, t), xx — Bx, t\ = Ct. 

where A < 1, B < 1 and C = B 2 ~ 2s A, so that the same equation will be satisfied 
after the change of scale. We propose here to do the same iteration for the solution 
u in terms of the variables x' and t. The equation will then take the modified form 



(12.4), that will be satisfied again by the iterates, just as it is written. It is true that 



the velocity v(t) will change from iteration to iteration according to the rule 

v 1 (t) = ^-v(Ct) = B 1 - 2s v(Ct), 
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which follows both from the geometric transformation, and from the definition of v 



in (12.2). Therefore, after the first geometrical transformation such repeated itera- 
tions conserve the same correspondence for all subsequent steps. In other words, the 
geometrical transformation done in the first step will hold for all remaining steps: if 
the set of coordinates at that moment is (x n ,t n ), we obtain a set of newly distorted 
coordinates (x' n ,t n ) by the formula 

x' n (t) = x n (t) - S n (t), S' v (t) = v n (t) 

then this is just a scaled version of the original transformation for n = 0. Summing 
up, since the contractions in the upper bound for u happen with a constant rate 
1 — \l in cylinders that shrink in space and time also with a fixed rate, we conclude 
in a standard way C a -regularity with respect to the transformed variables. But 
since the coordinate transformation is done only once and is Lipschitz continuous, 
this means the same type of Holder regularity for u with respect to the original 
coordinates (x,t). Of course, the Lipschitz constant of the transformation depends 
on Mi = sup t ||n(t)|| L i. 

The analysis of the second alternative is easier since we are converging along the 
iterations to an equation with constant diffusivity coefficient. We leave the easy 



details to the reader. This ends the proof of Theorem 5.1 for 1/2 < s < 2. □ 



13 Extension of the existence theory 

After these results, we can extend the existence theory to all nonnegative and inte- 
grable initial data. 

Theorem 13.1 For every uq G L 1 (M. n ), uq > 0, there exists a continuous weak 
solution of the FPME in the following sense: there exists a function u(x, t), continuous 
and nonnegative in Q = W n x (0,T) such that 

u E L°°(0, oo : L\R N ) n L°°(R N x (r, oo) for all r > , 

K{u) e L\0,T: W^(M. n )), uV)C(u) e L\Q T ) 

and the identity 

(13.1) u(rj t - V/C(w) • Vr?) dxdt + / u (x) rj(x, 0) dx = 
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holds for all continuously differentiable test functions rj in Qt that are compactly 
supported in the space variable and vanish near t = T. 

Proof, (i) We take a sequence of initial data uo n that are nonnegative, smooth and 
decaying at infinity as required in [II] . We assume that ||« ori ||i < C\ for all n > 1 
and uo n — > u in L 1 (IR Ar ). Then there exist solutions u n (x,t) and there are estimates 
like u„, > and 

u n (x,t) dx = \\u 0n \\i < Ci Vn > 1, t>0. 



ii) We also know from Theorem 4.1 that the solutions are bounded for t > r > 



with bound that depends only on the L 1 -norm of the initial data 

< u n (x,t) < C{n,s)\\u 0n \\Jt- a < CCjr a . 

We also know that uniformly bounded nonnegative solutions are C 13 smooth for some 
(3 > with uniform Holder constants, so that after passing to a subsequence we have 

u n (x, t) — > u(x, t) uniformly on compact subsets of Q . 

(iii) We now observe that for every t > we have u n (t) G L 1 R L°° which together 
with the Sobolev embedding (or Riesz embedding) gives 

H(u n {t)) e L r (R N ) V r > n/(n - s) 

and the embedding is compact into L^ oc (R N ). Indeed, if u n G L * p then 

1 _ 1 s 

r p n 

In practice we will need the estimate 



p> 1. 



H(u n (t)) 2 dx<C\\u n (t)\\l, 1 = --- 

F 2 p n 



But since 

\\u n (t)\\ p p < IKWIIilKWIlS- 1 < cwuml^-^r^, 

we get the decay estimate for this energy in the form 

H(u n (t)) 2 dx < C ||u(t) || J( 1 +Cp-i)'t)/J' t -2«Cp-i)/p = c 
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with 



2a(p - 1) 2n (n - 2s) n - 2s 

A — 1 — — 1 — -; : — 1 — 



p (n + 2-2s) 2n n + 2 - 2s n + 2 - 2s 

(iv) We also have for every ti > t\ > the energy inequality 

^ y H(u n {t 2 )) 2 dx + jf" y w n |Vi^ n | 2 cfcdt < i y H{u n (ti)fdx 

Therefore, we have a uniform estimate for the integral JJ u n \ VKu n \ 2 dxdt from t = 
r > up to t = infinity, and we get the same decay rate. 

(v) Next, we address the step of passing to the limit of the weak formulation of the 
solutions u n when the test function r\ G C X {Q) that compactly supported in Q away 
from t = 0: 

(13.2) J J u n {rj t — V/C(w n ) -Vrj) dxd + J Uo n (x) r)(x,0) dxt = 0. 

There is no problem in the convergence of J J u n rjt dxdt. If r\ is supported away from 
t = then the second term converges since we can use the compactness of the map 
u i — y VK(u) from L 1 R L°° into U for some r 

(vi) We now want to check that the initial data are taken. For that we show first that 
vSJKu G L\ t also near t = 0. Here is the calculation in a cylinder = M. N x (tk, tk-i) 
with 4 = 2~ k : 



[[ u\VK(u)\dxdt<(ff udxdt) 1/2 (ff u\V K(u)\ 2 dxdt) 

JjQ k JjQ k JjQ k 

u,Uk-i) 1,2 {\ J Hiu^ydx^KCt 1 ' 2 ^ 



< (IkJI,/,,-,) 1 "'' 2 ^ / HdKt,))^!.,-) 1 ^ < Ctl /2 t[ x - 1)/2 = Ct x k /2 . 

With this it is possible to pass to the limit in the term JJ u n V/C(n n ) ■ Vrjdxdt for rj 
that does not vanish near t = by estimating the integral for small t as uniformly 
small, and then proving the convergence for t > r > using the known regularity. 



This proves that the definition of solution according to Formula (13.1) is true. 

(vii) We can also check that the initial data are taken as traces. This depends on a 
uniform estimate of the difference u(t) — uq for small r in some adequate norm. We 
have that for a Cq test function rj(x) 



u{t) — uq\ rj dx < J J \d t u\r] dxdt < J u \VK(u)\ |Vt/| dxdt 
<C{jJ u\Vv\ 2 dxdt) 1/2 (JJ u\VK(u)\ 2 dxdt) 1/2 <Ct x/2 . 



52 



This is what makes u{t) — > uq in W loc ' in the limit. Or we may use J. Simon's 
compactness results [18], since we are actually proving that u t G L\{W~ 1,1 ). □ 
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